Deconfinement phase transition and the quark condensate 
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We study the dual quark condensate as a signal for the confinement-deconfinement phase tran- 
sition of QCD. This order parameter for center symmetry has been defined recently by Bilgici et 
al. within the framework of lattice QCD. In this work we determine the ordinary and the dual 
quark condensate with functional methods using a formulation of the Dyson-Schwinger equations 
for the quark propagator on a torus. The temperature dependence of these condensates serves to 
investigate the interplay between the chiral and deconfinement transitions of quenched QCD. 

PACS numbers: 



Introduction 

The chiral and deconfinement transition of QCD is a sub- 
ject of continuous interest. One of the unsolved problems 
is the question of an underlying mechanism relating both 
phenomena. Strictly speaking chiral and deconfinement 
phase transitions only occur in opposite sectors of the 
theory. The chiral condensate acts as order parameter 
for the chiral phase transition at vanishing quark masses, 
m = 0, whereas the Polyakov loop signals center symme- 
try breaking at the deconfinement transition for m oo. 
At intermediate masses both transitions develop a fasci- 
nating interplay [l[ , which is not yet understood in detail. 
Nonpcrturbative methods such as lattice gauge theory or 
functional methods are needed to explore this issue. 

A recent development that sheds light on this connec- 
tion is the investigation of spectral sums of the Dirac 
propagator and their behaviour under center transforma- 
tions. Initiated by a work of Gattringer spectral sums 
have been explored in Refs. [3, 0, S In particular it 
has been shown that the quark confinement mechanism 
is entirely encoded in the low lying spectral modes of the 
Dirac operator which are also responsible for chiral 
symmetry breaking via the celebrated Cashcr-Banks re- 
lation. The spectral sums constitute order parameters 
for the deconfinement transition of QCD 0, Q. 

In general, order parameters for deconfinement arc not 
easily accessible by functional methods. In Q a method 
has been developed to determine the Polyakov loop po- 
tential from the effective action, whereas in Q the ana- 
lytic structure of the quark propagator has been used to 
distinguish between the confined and dcconfincd phase. 
In this letter we report on a calculation of another order 
parameter for deconfinement, the dual quark condensate 
or 'dressed Polyakov loop' with functional methods. 

This dual condensate Ei is defined by the Fourier- 
transform 
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the temporal direction of Euclidean, quenched QCD. 
Note that the usual antiperiodic boundary conditions for 
fermions require if = tt, whereas ip = corresponds to 
periodic boundary conditions; here we vary tp in the inter- 
val [0, 27r]. In order to explain why the quantity Ei is of 
considerable interest we note that the (^-dependent quark 
condensate {'^p4')^p can be represented by a sum over all 
possible closed chains of link variables, i.e. closed loops 
I, on a lattice. One obtains: 
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of the ordinary quark condensate {'4>'4')ip evaluated us- 
ing J7(l)-valued boundary conditions with angle ip in 



where U (l) denotes the closed chains of links including 
some sign and normalisation factors, see Q for details. 
Each loop consists oi \l\ links and is weighted with the 
corresponding power of the inverse quark mass m. Each 
time such a closed loop winds around the temporal di- 
rection of the lattice it picks up a factor e^**^ from the 
?7(l)-valued boundary condition introduced above. This 
leads to a weighting e*'''"^'^ where n{l) is the winding 
number of a given loop I. The dressed Polyakov loop Si, 
Eq. ([1]), projects out loops with n{l) = 1. It transforms 
under center transformation in the same way as the con- 
ventional Polyakov loop [9| and is therefore an order pa- 
rameter for the deconfinement transition. The numerical 
agreement between dressed and conventional Polyakov 
loop has been estabhshed for gauge groups SU{3) [Toj 
and, remarkably, also for the centerless G(2) [ll|. 

In this work we determine the dual quark condensate 
from the Dyson-Schwinger equations (DSEs) of Landau 
gauge QCD at finite temperature • However, we wish 
to point out that our method is sufficiently general to 
be of equal use in other functional approaches as e.g. 
the functional renormalisation group [ij]. We evaluate 
the ordinary and the dual quark condensate from the 
trace of the quark propagator in a formulation of the 
DSEs on a torus. Our investigation of chiral symmetry 
breaking and deconfinement complements corresponding 
ones at vanishing temperature, see e.g. [l4| . and provides 
an interesting tool for further investigations of the QCD 
phase diagram by functional methods. 
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FIG. 1: The Dyson-Schwinger equation for the quark propa- 
gator. FiUed circles denote dressed propagators whereas the 
empty circle stands for the dressed quark-gluon vertex. 



Dyson-Schwinger equations on a torus 

We work in Euclidean space with compact time and space 
directions, i.e. box size V ~ L^x l/T with temperature T 
and 1/T << L. We choose periodic boimdary conditions 
in the three spatial directions for the quark and gluon 
fields. The gluon field also obeys periodic boundary con- 
ditions in the temporal direction. For the fcrmion field 



■0 we use the generalised, C/(l)-valued boundary condi- 
tion ip{x,l/T) = e^'^^{x,0) as described above. In loop 
integrals in momentum space this results in Matsubara 
modes ujp{nt) = {2'KT)(nt + </?/27r) in the p4-direction, 
which depend on the boundary angle ip € [0, 27r[. In 
the spatial directions we have the usual Matsubara sums 
which are treated with the techniques described in [l5l |. 

In Euclidean momentum space at nonzero temperature 
the renormalised dressed quark propagator is given by 

S{p,ujp) ^ [i-f4LL>pC{p,ujp) + i^A{p,ujp) + B{p,Up)]~^ , 

(3) 

with vector and scalar quark dressing functions C, A, B. 
For the bare quark propagator Sq we have A = C = 1 
and B ~ Zmm, with bare quark mass m. The DSE for 
the quark propagator, shown in Fig. [1] reads 
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where the sum is over temporal and spatial Matsubara 
modes. The Casimir Cp = {N^ — l)/-^c stems from the 
colour trace; in this work we only consider the gauge 
group SU{2). Furthermore, D^i, denotes the gluon prop- 
agator in Landau gauge and F^ the (reduced) quark- 
gluon vertex. The renormalisation factors Z2, Z„i and 
Z3 are determined in the renormalisation process. 

In order to solve this equation we have to specify ex- 
plicit expressions for the gluon propagator and the quark- 
gluon vertex. For the momentum range relevant for 
Eq. ([4]) we nowadays have very accurate solutions for the 
gluon propagator at zero temperature from both, lattice 
calculations and functional methods, see e.g. \m and 
references therein. The temperature dependence of the 
gluon propagator, however, is much less explored. In 
|17| a combined lattice and DSE study records a different 
temperature dependence of the electric and magnetic sec- 
tor. Whereas the magnetic part of the propagator seems 
to be indifferent to the deconfinement phase transition, 
the electric part is strongly increased at and around the 
critical temperature « 300 MeV. Although the lat- 
tice data still have considerable systematic errors they 
may correctly represent the qualitative temperature de- 
pendence of the gluon. We therefore use a temperature 
dependent fit to the data given by 



T[MeV] 





119 


298 


597 


aT(r) 


1 


1 


1.34 


1.65 




1 


1 


0.8 


4.0 



TABLE I: Temperature dependent fit parameter of Eq.(I7}. 
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with q = {q,^q) and dressing functions ZT{q,(jJq) and 
ZL^q^LOq). The transverse and longitudinal projectors 
with respect to the heat bath are given by 

Pi:Aq) - (6) 

with (i, j = 1 . . . 3). The SU{2) lattice resuhs of ReL [13 
arc well fitted by 



ZT,L{q,u}q,T) = 
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with the temperature independent scale A = 1.4 GeV and 
the coefficient c = 9.8 GeV^. For gauge group SU{2) we 
have /3o = 22/3 and 7 = —13/22 in the quenched theory 
and we renormalise at a(fi) ~ 0.3. The temperature de- 
pendent scale modification parameters aT.L{T) are given 
in table m In order to extend this fit to temperatures not 
given in the tabic we assume aT,L{T) to be temperature 
independent below T = 119 MeV and only slowly rising 
above T = 597 MeV. For T e [119,597] MeV we use 
cubic splines to interpolate smoothly between the val- 
ues given in table HI We expect the systematic error of 
this procedure to be of the same order as the systematic 
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errors inherent in the lattice data. We also inherit the For the quark-gluon vertex with gluon momentum q = 
scale determined on the lattice using the string tension {q,ujq) and the quark momenta p = {p,ujp),k = {k,u!k) 
= 0.44 GeV we employ the following temperature dependent model 



r,,(q,fc,p) = Zg ( (54^74 \-Oj^Jj 
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where 5 = —9/44 is the anomalous dimension of the ver- 
tex. Note that because of 7 + 2(5 = — 1 the gluon dress- 
ing function together with the quark-gluon vertex behave 
like a running coupling at large momenta; this is a neces- 
sary boundary condition for any model interaction in the 
quark DSE. The dependence of the vertex on the quark 
dressing functions A and C is dictated by the Slavnov- 
Taylor identity. The remaining fit function is purely phe- 
nomenological, see e.g. where an elaborate version of 
such an ansatz has been used to describe meson observ- 
ables. Here we use di = 7.6 GeV^ and da = 0.5 GeV^ A 
variation of these parameters shifts the critical tempera- 
tures of both, the chiral and the deconfinement transition 
but leaves all qualitative aspects of the results presented 
below unchanged. This is also true if one changes the 
deep infrared behaviour of the vertex to the (infrared 
divergent) form extracted in [l^ and used in the 
confinement-deconfinement phase transition is insensitive 
to the question of scaling vs. decoupling in the sense 
specified in [l^. Details will be reported elsewhere. 

The DSE is solved numerically using rit = 8 Matsubara 
frequencies and a three- volume of V = (Sfm)"^, which is 
large enough to avoid any significant volume effects. The 
renormalisation conditions are C(^) = 1 and B{p) = m 
with /J, = (/I, ttT) and = 20 GeV^. The resulting quark 
dressing functions are subsequently used to determine the 
quark condensate according to 



AZ2NcT 
~L3 



E 



52 



(9) 



where we indicated the dependence of the frequencies on 
the generalised ?7(l)-boundary conditions. For nonvan- 
ishing bare quark masses m this expression is divergent 
and, at least in the continuum limit, has to be renor- 
malised accordingly. For the purpose of this letter, how- 
ever, is is sufficient to work with the regularised expres- 
sion at fixed ultraviolet cut-off corresponding to fixed lat- 
tice spacing. 

Numerical results 
First we explore the dependence of the quark condensate 
on the boundary angle (p which can be read off Eq. 
Each loop winding n times around the temporal direction 
of the torus contributes a factor cos{nip). Consequently 
the integrand in Eq. ([1]) can be expanded as a series in 
cos(niy9) and is symmetric in the interval [0, 27r]. This 



behaviour is clearly seen also in our numerical result for 
= (^ip)^ shown in FigH For T = 200 MeV far below 
the deconfinement transition we find almost no angular 
dependence of the condensate. This is especially true 
for the heavier quark mass m — 60McV. For T = 400 
MeV far above the transition the behaviour is markedly 
different and we observe the characteristic cosine-type of 
behaviour of the condensate that we expect from Eq. ([2]). 
This is in nice agreement with the results of Ref. [a] on 
the lattice. The angular dependence of the condensate 
can be fitted well with a series in cos{n(fi). Whereas for 
the heavier mass m = 60 MeV terms with n < 3 are 
sufRcient one needs at least terms up to n < 7 for the 
smaller mass m = 10 MeV. Thus the closer one gets 
to the chiral limit the more contributions from Polyakov 
loops with higher winding number around the temporal 
direction are significant. This is a direct consequence of 
the mass factor l/ml'l in the expansion Eq. ^ and also 
seen on the lattice [l^ . 

By far the largest contribution to the angular depen- 
dence of the condensate comes from loops with n = 1 
which are projected out by the Fourier transformation 
(H]) to the dual condensate. The resulting temperature 
dependence of the dressed Polyakov loop is shown in 
Fig. [3] together with the conventional quark condensate 
and their derivatives. One clearly observes a change of 
behaviour in both, the conventional and the dual condcn- 
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FIG. 2: The angular dependence of the quark condensate 
= (V'V')*' below and above the deconfinement transition 
for two different quark masses. 
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FIG. 3: The temperature dependence of the dressed Polyakov- 
loop El and the conventional quark condensate = 
{ipij)) ip=-n together with their derivatives for m — 10 MeV. 

sate above T = 270 MeV. The temperature derivative of 
both quantities has a peak in the region Tc = 300 — 320 
MeV signaUing the chiral and deconfinement transition. 
Note that although both transitions are calculated from 
quantities with direct relation to the spectral properties 
of the quark propagator they do not necessarily give the 
same transition temperatures. On the contrary, the chi- 
ral transition occurs about 10 — 20 MeV below the decon- 
finement transition. Whether the quantitative aspects of 
this difference is a model-independent result has to be 
investigated in more detail. 

Finally we wish to point out that the absolute magni- 
tude of our dressed Polyakov loop is about a factor of five 
smaller than the one calculated on the lattice Since 
the condensate and subsequently its dual is a renormali- 
sation point dependent quantity we attribute this differ- 



ence to a yet undetermined renormalisation factor. This 
issue has to be clarified in future work. Furthermore we 
note that the dressed Polykov loop is not strictly zero 
in the low temperature region as expected in quenched 
QCD. This is certainly a consequence of our vertex ansatz 
which does not strictly represent quenched QCD but in- 
cludes a (small) amount of unquenching effects. 

Concluding remarks 

In this letter we have determined the dual condensate or 
dressed Polyakov loop by solving the Dyson-Schwinger 
equations for the quark propagator on a compact mani- 
fold. This order parameter for center symmetry measures 
contributions from closed loops winding once around 
the time direction. We observe a significant rise of the 
dressed Polyakov loop around and above the deconfine- 
ment transition temperature together with a significant 
decrease of the ordinary quark condensate. The angu- 
lar dependence of the quark condensate shows a char- 
acteristic dependence of the quark mass: the lighter the 
quark the more contributions arise from loops with wind- 
ing numbers larger than one. An obvious next step is to 
investigate what happens to the transition temperatures 
when the backreaction of the quarks onto the Yang-Mills 
sector is included. Particularly interesting in this respect 
is the case Nf =2-1-1 which recently is a subject of 
intense debate 0, [2l| in the lattice community. 

Acknowledgments 

I thank Falk Bruckmann, Christoph Gattringer, Jens 
Miiller and Jan Pawlowski for discussions. I am grateful 
to Axel Maas for discussions and for making the lattice 
data of Ref. [13 available. This work has been supported 
by the Helmholtz-University Young Investigator Grant 
number VH-NG-332. 



[1] F. Karsch, 'arXiv:hep-lat /9903031 ' 

[2] C. Gattringer, Phys. Rev. Lett. 97 (2006) 032003. 

[3] F. Bruckmann, C. Gattringer and C. Hagen, Phys. Lett. 

B 647 (2007) 56. 
[4] F. Synatschke, A. Wipf and C. Wozar, Phys. Rev. D 75 

(2007) 114003; 

[5] E. Bilgici, F. Bruckmann, C. Gattringer and C. Hagen, 

Phys. Rev. D 77 (2008) 094007. 
[6] F. Synatschke, A. Wipf and K. Langfeld, Phys. Rev. D 

77 (2008) 114018. 

[7] J. Braun, H. Gies and J. M. Pawlowski, 'arXiv:0708.2413' 

[hep-th]; F. Marhauser and J. M. Pawlowski, 
'arXiv:08r2j[l44 [hep-ph]. 
[8] A. Bender, D. Blaschke, Y. Kalinovsky and 

C. D. Roberts, Phys. Rev. Lett. 77 (1996) 3724. 
[9] A. M. Polyakov, Phys. Lett. B 72 (1978) 477; 

L. Susskind, Phys. Rev. D 20, 2610 (1979). 
[10] F. Bruckmann, C. Hagen, E. Bilgici and C. Gattringer, 

PoS LATTICE2008 (2008) 262. 
[11] J. Danzer, C. Gattringer and A. Maas, JHEP 0901 

(2009) 024. 



[12] C. D. Roberts and S. M. Schmidt, Prog. Part. Nucl. Phys. 

45, SI (2000); A. Maas, J. Wambach and R. Alkofer, Eur. 

Phys. J. C 42 (2005) 93. 
[13] J. Braun and H. Gies, JHEP 0606 (2006) 024. 
[14] R. Alkofer, C. S. Fischer, F. J. Llanes-Estrada and 

K. Schwenzer, Annals Phys. 324 (2009) 106. 
[15] C. S. Fischer, R. Alkofer and H. Reinhardt, Phys. Rev. D 

65 (2002) 094008; C. S. Fischer and M. R. Pennington, 

Phys. Rev. D 73 (2006) 034029; G. S. Fischer, A. Maas, 

J. M. Pawlowski and L. von Smekal, Annals Phys. 322 

(2007) 2916. 

[16] C. S. Fischer, A. Maas and J. M. Pawlowski, 

arXiv:0810.1987 [hep-ph]. 
[17] A. Cucchieri, A. Maas and T. Mendes, Phys. Rev. D 75, 

076003 (2007). 

[18] G. S. Fischer and R. WiUiams, Phys. Rev. D 78 (2008) 
074006. 

[19] F. Bruckmann, G. Gattringer and G. Hagen, private com- 
munication. 
[20] Y. Aoki et ai, arXiv:0903.4155 [hep-lat]. 
[21] A. Bazavov et o/.. .arXiv:0903.4379. [hep-lat]. 



